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Based on a tight-binding approximation, we present analytical solutions for the wavefunction 
and propagation velocity of an electron in armchair graphene ribbons. The derived expressions 
are used for computing the transmission coefficients through step-like and barrier-like potentials. 
Our analytical solutions predict a new kind of transmission resonances for one-mode propagation 
in semiconducting ribbons. Contrary to the Klein paradox in graphene, this approach shows that 
backscattering for gapless mode is possible. In consistence with a higher order k • p method, the 
backscattering probabilities vary with the square of the applied potential in the low-energy limit. 
We also demonstrate that gapless-mode propagation through a potential step in armchair ribbons 
can be described by the same through-step relation as that for an undimerized ID chain of identical 
atoms. 
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I. INTRODUCTION 

Recently graphene sheet, a monolayer of covalent-bond carbon atoms forming a dense honeycomb crystal, has been 
obtained experimentally^"^. Even at room temperature submicron graphene structures act as a high- mobility electron 
or hole conductors^. This phenomenon has brought to life a promising field of carbon-based nanoclcctronics, where 
graphene ribbons (GRs) could be used as connections in nanodevices. 

The presence of armchair- and zigzag-shaped edges in graphene has strong implications for the spectrum of tt- 
electrons^i^i^ii and drastically changes the conducting properties of GRs. Especially, a zigzag edge provides localized 
edge states close to the Fermi level {Ep = 0), which leads to the metallic type of conductivity. In contrast, any 
localized state does not appear in an armchair GR. An armchair GR can be easily made to be either metallic or 
semiconducting by controlling the width of the current channel. 

Most of the intriguing electronic properties of GRs were first predicted by the tight-binding model^'^'^'^'^ i^^'^^ . 
These features are also well reproduced in the k • p methodi^iii^ based on the decomposition of linear Schrodinger 
equations in the vicinity of zero-energy point. The commonly used k • p equation is a two-dimensional analog of the 
relativistic Dirac equation which is a certain approximation of general tight-binding Schrodinger equations. 
The continuum Dirac description of the electronic states near E — Ep has been shown to be quite accurate by 
comparing with a numerical solution in the nearest-neighbor tight-binding model^i^. Particularly, in the framework 
of the Dirac approach, perfect electron transmission is predicted through any high and wide potential barriers, known 
as the Klein paradox in graphenai^. 

Recent analytical studies based on the relativistic Dirac equation have considered transmission through barriers^, 
graphene wells^iii, double-barriersi^, quantum dotsi^, superlattices^S, and n — p junctions^!. However, transport 
properties of these structures in the tight-binding model has not been analyzed systematically yet. Considerable 
attention has been paid to the effects of conductance quantization in GRs in the absence of scattering, where the 
quantization steps clearly indicate the number of propagating states crossing the Fermi leveUi. 

In this paper, we develop an exact model of electron transport through armchair GRs based on the nearest-neighbor 
tight-binding approximation. (In contrast to the zigzag GRs, the armchair GR spectrum does not have a special band 
of edge states complicating the description.) We obtain simple analytical expressions for normalized wavefunctions 
and propagation velocity of an electron wave in armchair ribbons. No similar explicit solutions have been obtained for 
armchair GRs yet. In a tight-binding study by Zheng^ of the electronic structure in armchair ribbons, the derived 
analytical form of wavefunction has not been used for describing charge transport in GRs. 

On the basis of the analytical solutions for the wavefunction and propagation velocity of an electron in an infinite 
armchair GR, we describe the scattering problem through step-like and barrier-like profiles of site energy along the 
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ribbon in closed form. Solving the corresponding equations for the scattering amplitudes, we obtain exact transmission 
coefficients in the atomistic (tight-binding) model analytically. For the energy region close to the Fermi level we 
compare our analytical expressions with their analogs obtained in the continuous (Dirac) model^^. Our general 
conclusion is that the validity of the relativistic approach based on the standard k • p equation^'^ is overestimated for 
the description of charge transport in armchair GRs. Specifically, for one-mode propagation in semimetal GRs, the 
deviation of transmission coefficients T from the unity near the zero-energy point is proportional to the square of the 
step (barrier) magnitude, which contradicts the known Klein paradox in graphenel^,. It is worth noting that similar 
deviations of transmission coefficients are intrinsic for the extended k • p equation, based on higher order k • p terms 
in the low-energy limit. Such high-order approximations lead to predicting small backscattering due to the effects of 
trigonal warpingi^ of the bands, destroying perfect transmission in the channel. 

The paper is organized as follows. In the framework of the nearest-neighbor tight-binding formalism we calculate 
energy dispersion (SecHJ) and eigenstates (Sec. IIIip of electrons in the honeycomb graphene lattice. The lattice, in 
contrast to Ref;22,, is considered as a set of rectangular elementary cells of four nonequivalent atoms (Fig. [T]). This 
assumption simplifies the description of electron transport processes in GRs in the presence of step-like or barrier-like 
site-energy profiles. Following Ref.^^, in Section llVI we derive boundary conditions for infinite GRs with armchair 
edges. By using a combination of early found eigenstates of the honeycomb lattice, we obtain normalized wave 
solutions for electrons in armchair GRs and formulate the quantization law for the transverse component fc™ of the 
wave vector. As a result, the full 7r-electron band spectra of graphene ribbon breaks into a set of subbands, whose 
features are discussed in Section IVl 

In Section IVTl we find the group velocity and propagation direction of an electron in an armchair ribbon. It is shown 
that the group velocity is proportional to sin 0, where 6 is the electron phase shift between two neighboring unit cells 
of the ribbon. In sections IVIII and I Villi the transmission probabilities through step- like and barrier-like potentials are 
calculated. Based on the exact solutions, in Section [IXl we demonstrate that the transmission coefficients exhibit the 
new kind of resonance occurring when cos9 = cos 9, where 9 is the inter-cell electron phase shift in the region where 
the potential exists. The number and positions of the transmission resonances are shown to be different in the regions 
fc™ < 2tt/3 and fc™ > 2tt/3. As for one- mode propagation in semimetal GRs (fc™ = 27r/3), our theory predicts no 
unit propagation. In Section |Xl by expanding the expressions for cos 9, cos 9 in the vicinity of -E = 0, we explain why 
the Klein paradox does not hold in armchair GRs. 

In Sec. IXI[ we discuss through-step and through-barrier probabilities and provide some examples. Conclusions are 
presented in Sec. IXIII Appendices [Xl f on electron flux in the tight-binding model) and[Bl(on through-step transmission 
in a linear chain of identical atoms) describe some theoretical tools used in the paper. Appendix [C] contains some 
technical details related to obtaining our approximate solution from the Dirac solutionis. 



II. DISPERSION RELATION FOR HONEYCOMB GRAPHENE LATTICE 

We consider the graphene honeycomb lattice as a set of rectangular elementary cells of four atoms a = I, X, p,r (see 
Figure [T]). By taking the tight-binding representation for molecular orbitals 

oo 

n,m— — oo a 

we come to the set of linear Schrodinger equations 

'—E'4>n,m,l = i'n,m,X + V'n,m+1,A + ^n-l,m,r, 
-Elpn,m,\ 
-E')pn,m,p 

— '0n,m,p ~l~ '4^n,'m-\- l,p + ^n+l,m,l, 

with respect to wave function components ipn,m.a = (^|?t., m, a). Here |n, m, a) is 2pz orbital of a-th atom in 
{n, m}-elementary cell, E = E/\(3\ is the electron energy in units of |/3|, /3 < is a transfer integral between the 
nearest-neighbor carbon atoms. Site-energies of carbons equal zero and serve as a reference. 
We look for a solution of Eq. ([T]) in the form 

V'n,m,a = 0ae^'""+*'='"'", a = l,X,p,r, (2) 

where the dimensionless wave numbers k"{= 3kxac-c), fc'"(= \/3kyac-c), and co-factors (j)a are unknown coefficients, 
kx and ky are the components of the wave vector along n and m directions respectively. The lengths 3ac-c and 
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FIG. 1: Sketch of 2D honeycomb lattice. The central shadowed block corresponds to an elementary cell {n, m} consisting of 4 
carbon atoms, labeled as I, A, p, r. The lengths 3ac-c and \/3ac-c are the translation periods in n and m directions, ac-c is 
C - C bond length. 

-\/3ac-c are the translation periods of the graphene sheet as depicted on the Fig. [l] where ac-c is the C-C bond 
length. 

Plugging in the solutions ^ into the system ^ gives the system of linear equations 

'E(t)i + (1 + e''^") + e-*'=^'>,. = 0, 
(l + e-*'")0j + S0A + 0p = O, 

^A + £;0p + (l + e-»^'")0, = 0, ^ > 

^e*'=>, + (l + e'^'")0p + ^0, = 0. 



A nontrivial solution of Eq. ([3]) demands a zero determinant. It leads to the dispersion relation 

(4) 



r£;2 _ 1 _ 4cos2 —1^ 

cos fc" = /(£;, fc™), /(i?,fc-) = l 2 J 

8 cos^ ^ 



defining the non-dimensional wave number /c" in terms of k"^ and electron energy E. It follows from Eq.Q that 

i?^ = l±4cos — cos |-4cos^ — , (5) 

2 2 2 ^ ' 

where index ± distinguishes two branches of the dispersion relation. 

III. WAVE SOLUTION IN GRAPHENE LATTICE 

If the relation ^ holds, system ^ becomes degenerated, and we can express the coefficients 0/, <j>\, (pp in terms 
of (jjr, omitting the fourth equation in ([3]). The result can be written as follows 

(j>a^C<fa, a^l,X,p,r, (6) 

where 

(£;2_i)e-*'=" +4cos2 4^ ~ + (1 + 6-''=") 

A. - \ ' ^ -tlx - ^ ' ^ ' 



£;(£;2-l-4cos2^) ' ^ £;2_l_4cos2^ 



i;(i;2_i_4cos2 ^) 
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FIG. 2: Sketch of infinite graphene ribbon with armchair edges. The soUd-framed block is the unit cell n of GR. 



Exploiting the energy dispersion relation (|4|), one can see that — 1 for the real- valued fc". Introducing a new 
function 

.„ £;2 - 1 + 4e*'=" cos2 ^ l±2cos4^e*'="/2 

e = — — f7) 

S(S2_i_4cos2*^) E ' ^ ' 

we can get 

0, = e'(''-'="), = ±6-^^'='"+"")/^ 0^ = Te-K'='"+'=")/2+^«, ^^^^^ (8) 

The choice of upper/lower signs in (O, dH) is determined by the sign in ([5]). Plugging in Eq. ([8]) into Eq. ([6]), and 
then the resulting expression into ^ gives the expressions for electron eigenstates in the honeycomb lattice 
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_^gifc"(«-l/2)+ifc"(m-l/2)^ a = A, 

-f-c; , ct 



IV. WAVE SOLUTION FOR ARMCHAIR GRS 



Unlike the 2D honeycomb lattice, the electron wave function components Tpn,m,a in GRs with armchair edges (see 
Figure 12) are the solution to the linear Schrodinger equations ([T]) only for inner elementary cells (1 < to < M). For 
boundary cells with m= \ ov m = M + 1, one needs to take into account that boundary carbons have only two 
neighbors (see Fig. [2]). It demands the wave function to vanish on the set of absent sites^ i^^i^^ 



Similar boundary conditions were used in Refs^i^S.. 

To satisfy the relations pO|) . we represent the solution as a linear combination of the states ([9]), 



V'n,m,a(^ ) — ^n,m,a(^ !^ ) '0n,m,a(^ 

Since the definition ([7]) implies that 6'(— fc™) — d{k™), we finally obtain 



'4'n,m,a{k ) — C < 



(^tk"{n-i)+te sinfc™TO, 

_^g»/c"(n-l/2) sinfc™(TO- 1/2), 



e*'^"" sinfc™m, a = 



a ^ I, 
a — X, 



(10) 



(11) 
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where 

■Kj 



j = l,2,...,AA, (12) 

is the set of discretized transversal wave numbers, and the unknown constant C"(= 2iC) can be found from the 
normahzation condition over a unit ceU of the armchair GR, 



m— 1 m— 1 



m— 1 

Then we finaUy get C = [2(A/' + l)]-i/2 



V. BAND SPECTRA OF ARMCHAIR GRS 



Due to the transverse momentum quantization (|12p . the fuU 7r-electron band spectra of an armchair GR breaks 
into a set of subbands connected with each mode j independently. Given E and fc™ = fc™, the dispersion relation 
(|4|) allows us to establish some important peculiarities of the j-th part of the band spectra. Following the theory of 
a one-dimensional crystal-'* with an arbitrary electronic structure of the elementary cell, each j-ih band of graphene 
electron spectra includes 4 subbands, symmetrically disposed with respect to i? = 0. Real values of longitudinal 
wave number fc" e [0, tt] define the region of propagating electron states (when [/(E', fc™)| < 1 in Eq.Q). The cases 
/(i?, fe™) > 1 and /(i?, fc™) < — 1 are related to the complex values of fc™ = i5 and fc™ = tt + ii5, respectively. They 
refer to the forbidden zones - the gaps between neighboring electron subbands and the regions above the highest 
subband and below the lowest one. The subband boundaries correspond to fc™ = or fc™ = tt and satisfy the 
solutions to f{E,k™') — 1 or f{E,k'"'') — —1, respectively. Since Eq.(l4|) defines only one fc™ for each energy level E, 
the subbands of the jth band cannot overlap, however, they may touch each other along the frontiers with the same 
fc'"(= 0,7r). 

The results of the band spectrum modeling for armchair GRs with J\f = 10—12 are represented in Figure[3l Similar 
shapes are observed for GRs with arbitrary number J\f. The common features of the band spectra are: (I) bands 
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FIG. 3: Examples of energy bands for armchair graphene ribbons with A/" = 10 — 12. The transverse momentum quantization 
results in subband series connected with each mode j independently. Gray (black) rectangles mark subbands, whose bottom 
corresponds to fc" = {k" = tt) and whose top - to fc" = tt (fc" = 0), as depicted in the framed panel. The border between 
black and gray rectangles in the positive- and negative-energy spectra determines the subband connecting fc" = vr frontiers, this 
is a degenerated solution to f{E,kJ^) = —1. The fundamental gap of a graphene ribbon is the minimal gap value A between 
the valence and conduction bands among all modes j — 1, . . . , A/". 



are bounded by the energy interval \E\ < 3; (II) the band structure is symmetric with respect to i? = 0; (III) if the 
number A/'-f 1 is divisible by 3, the mode j = 2{Af + l)/3 (or fc™ = 27r/3) does not possess the energetic gap between 
positive- and negative-energy bands, and, therefore, it possesses the semimetal type of electron conductivity-ii'^'2^. 
For other j-modes there is a gap 



A 



1 — 2 cos ■ 



(13) 



between positive- and negative-energy subbands, which corresponds to a semiconducting ribbon whose gap decreases 
with increasing J\f. 
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VI. PROPAGATION VELOCITY 

To obtain propagation velocity of the wave, we plug in the solution (fTTj) into the relation (A2) for the electron flux 
J. Then J — sin6'/2?i. Using the definition ([7]) and the first derivative of the dispersion relation ^ with respect 
to E, one can get 

j_ \f3\ 2sin/c"cos2^ _ \l3\ dE 

Representing the flux J as where w(/c") is the dimensionless group velocity, we flnally obtain 

, , „ , dE sin 9 
^; fc" = — = — — . 14 
^ ^ dfc" 2 ^ ' 

The sign of v in Eq. (|14p determines the moving direction of the propagating wave (|11|) . For a subband (=5), 
whose bottom corresponds to /c" = and top to fc" = tt (see Fig. |3]), the wave number fc" increases along with i? 
[dE/dk"^ > 0). Therefore, w > and the state pT|) describes the right-moving mode. 

Similarly, in a subband of type we have dE/dk" < 0, and the wave solution pT|) propagates to the left. 

Replacing /c" with — fc" in ([7]), pT|). we obtain the expression describing a right-moving mode in these subbands 

{g_,fc"(„_i)„,e sinfc^m, a = /, 

±e-*'"<""'/'^ sinfc™(m-l/2), a = A, 

^g-*/c"(„-i/2)-»9 sinfc™(TO - 1/2), a^p, ^^^> 

^-ik"n ginfc™m, a — r. 

From pT|) and (fT5|) . one can understand the physical meaning of the quantity 6. Expressing ipn+i,m,i{ik") via 
V'n,m,r(±fc"), wc get ipn+i,m,i{ik'^) = e^*^Vn,m,r(±fc"), i.e., an clcctron, propagating through the corresponding 
subbands, gains phase shift ±6* between two neighboring unit cells of GR. The phase shift plays an important role 
in electron transmission through the ribbons, as discussed below. Also note that sin 6* > inside the (=5)-subbands 
and sin 9 < in the subbands {—'i), which is also seen from the relation (fT4)l . 

The derived analytical form of the wavefunctions PT|) . and the expression ((TJ]) for propagation velocity will be 
used below to investigate transport properties of armchair GRs in the presence of step-like or barrier-like potentials. 

VII. TRANSMISSION PROBABILITY THROUGH A POTENTIAL STEP 

In this section, we derive the transmission coefficient for an electron with energy E, transverse and longitudinal 
wave numbers fc™, fc", meeting the electrostatic potential 

[Vo, n>l, 

in an armchair ribbon. We denote the wavefunction by ■(/)„ „ ^ for n < and V'n'm a ^'^^ n > 1. Let Uq = eVo/|/3|. 

The Schrodinger equations for all atoms in both leads are used to obtain the probabilities of electron scattering 
through the potential step. The comparison of these equations near the interface between two regions {n = 0, 1) 
allows us to write down the matching conditions'^ as follows 

left right left right 

V'o.m.r = V'0,m,r: "01,™,; = V"!,™,; , m=l,...,7V. (16) 

Similar conditions for the tight-binding model were used in Ref.^*^ and for the Dirac equation in Refsj ^^i^^'^^ . 

Since the matching conditions (|16p do not mix the modes, we need to sew only the j-th solutions in each region 
under study. Wave solutions, describing incident and reflected electron waves in the left lead, can be constructed from 
the states (dJ) and For the right region, one needs to take into account the shift of site energies produced by 

the applied potential. Introducing additional notation for longitudinal wave number /c", phase and group velocity 
v{k") in the right lead, 

fc" = fc"(i;^£;-?7o), e = e{E ^ E-Uo^k"" ^k""), w(fc") = ^, (i?) 
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FIG. 4; Sketch of a graphene ribbon with armchair edges modified by a step-like potential. Unit cells with n < and n > 1 
belong to the left and right leads of GR, respectively. 
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FIG. 5: Three possible types of electron subband disposition in the right region. The Fermi level is assumed to be located 
in a (=o)"Subband of the left lead and incident wave is described by the solution a) Electron propagates through a 

(=o)-subband in the right region; b) the Fermi level belongs to the band-gap of the right lead, which results in total reflection; 
c) electron propagates through a (^=5^) -subband in the right region. This case corresponds to the solution (|15|l for an outgoing 
wave in the lead. 



we first consider the case of propagation through the (=5)-subbands in both regions, as depicted on the a-panel of 
Fig. [5l The solutions describing the incident and reflected electron waves in the left lead and the outgoing wave in 
the right region should be taken in the foUowing form 



left 

n,rri.l 



sinfc™m, ip^ 



+ re 



sin k"^m, 



(18) 



Vn,m,l 



sin fcf TO, CI'mV = ^e'^"" sin fc™m, 



(19) 



where r and t are the amphtudes of the reflected and transmitted waves excited by the incident wave with transverse 
mode j. (We omit the factor C" for simplicity.) After plugging in [T5)) . ITOl) into the conditions we get the 

following two equations 

l + r = t, 



Their solution is 



t 



2i sint 



(20) 



The reflection and transmission probabilities R and T can be directly found from the scattering amplitudes r, t, 
using their definitions: R = |rp and T = |ipu(fc")/u(fc") = jipsin^/ sin 6*. In this way, we obtain the expression for 
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an electron propagating via (=J)-subbands in both leads 

A sin 9 sm0 



T = 



4 sin 9 sin 9 + (sin 9 — sin 9^ + (cos 9 — cos 9) ^ 



The b-panel of Figure [5] corresponds to the case when the electron energy belongs to the forbidden band in the 
right region (it is impossible for one- mode propagation through a semimetal ribbon with a zero-gap, fc™ = 27r/3). In 
this case, the wave number fc" takes imaginary values, and the expression for t in (|20p is to be omitted. At the same 
time, for imaginary-valued the expression for e*^ becomes real, which results in i? = 1 after using the r-solution 
from (pni). 

For an electron moving through a (=^)-subband in the right lead, the c-panel of Figure [5l we use the solution (jl5|) 
for an outgoing wave (see the definitions (fT7| ). Thus, we obtain 

—4 sin 6* sin^ 

T = 



—4 sin 9 sin 9 + (sin 9 -\- sin 6*) -I- (cos 9 — cos ( 



where sin0 < 0, as stated in Sec lVIl In a similar way, one can determine the transmission coefficient for an electron 
propagating through arbitrary subbands in the left and right leads 

ystop ^ 4 1 sing sin^l ^^^^ 

4| sing sing| -I- (| sing| - | sin 6*1)^ + (cos g - cos g)^ 

The common relation (|2ip must be supplemented with the expressions for cosg, cos 9 determined from ([7]) and (|17p. 
Then, from the energy dispersion ([¥]), we finally obtain 

cos9 ^ g{E,k"'), cos9 ^ g{E ~Uo,k"^), (22) 

where the function g is described by the simple relation 

+ 1 - 4cos2 

g{E, = — ^. (23) 

The analytical expression (PT|) with the definitions (P^ . is one of the central points of this paper since it makes 
it possible to compute the transmission coefficient through a potential step without recourse to the initial dispersion 
relation (0]). 

It is worth noting how similar the solution (PTjl to the relation (B4) is, though the relation (B4) is for through-step 
transmission coefficient in the ID undimerized chain of identical atoms. In other words, the inter-cell variables g, 9 
for GRs act as wave numbers fc, k in the linear chain model, whose unit cell degenerates to a zero-sized atom. We'll 
return to this issue in Sec. IIXI 



VIII. TRANSMISSION COEFFICIENT THROUGH A POTENTIAL BARRIER 

To obtain the exact expression for the electron transmission probability through a potential barrier 

n < 0, 

V{n) ^ {Vq, l<n<N, 
n > iV + 1, 

(see Fig. IH]), we denote the wavefunction by V'^^'^n „ for the left lead (n < 0), '4'nm a for the "in" region (1 < n < N) 

perturbed by the potential, V'n'^m a for the right lead {n > A'^ + l), and keep the same notation for the disturbed region, 
as in the previous section. By analogy with the matching conditions (|16p . we need to have the following relations 
satisfied 

in _ riglrt in _ riglit [^'i) 



As before, the solutions are determined using the states ([TT|) and 
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FIG. 6: Sketch of a graphene ribbon with armchair edges modified by the barrier-hke potentiaL W and D are the ribbon width 
and length of the disturbed region, respectively. 



For an electron propagating through a {—q) subband in the left lead the wavefunction expressions are presented by 
(fT5)) . The outgoing wave in the right lead is constructed from the state pT|) . 



n,m,a — tsm{kj m) <^ ^ik"n a = r 



ik" {n~l)+i9 ^g-ifc"(n-l)-i6l^ ^ _ ^ 



J,right _ 
r n,m, 

For the in-region we determine a solution 

! ae 

'^n,m,a = sin(fcj"m 

where a and 6 are additional unknown constants. The solution does not depend on the subband type of propagation 
in the disturbed region. The case of tunneling through a gap, with fc" imaginary, is not considered in the paper. 
Using the matching conditions (j24p . wc obtain the following four equations with four unknowns r, i, a and h 

ae»fc"JV ^ lJ^-^k-N^ (25) 





- re * 




'AT _ 







Solving with respect to i, we obtain 
where the functions 



sinA;"A^ sin^ 
" sin(fc"7V + 9) ' ^"'^ " sin(fc"Af + 9) ^ ' 

depend only on the "disturbed" variables fc" and 9. Since, by definition, T — the through-barrier transmission 
probability can be written in the following compact form^^ 

^2 n ^,-2 7) 



sin^e'sin^^+ (cos6' - cos^)^ sin^(fc"A^) 

It is easy to verify that the transmission coefficient (P7)) remains unchanged for propagation through a (=5^) -subband 
in external (left and right) leads. 

As follows from Eq. ^7} . the unit transmission occurs under the coincidence of cos0 and cos^. This can be viewed 
as a new type of resonance, which differs from the familiar resonance condition sin(fc" A'^) — 0. Some of its properties 
are studied in the next section. 
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FIG. 7: Behaviour of g{E, fc™) inside the positive- and negative-energy bands for A;™ = 2-k/3 (sohd line) and A:™ — 27r/3 =F 0.1. 
The shaded rectangles represent the band positions for the transverse wave number fc"" = 2n/3 — 0.1, where A is the energy 
gap for given fc™. 



IX. UNIT PROPAGATION CONDITION 



The transmission coefficients (pTj) and ((27|) through step-like and barrier-like potentials exhibit the unit propagation 
at 9s such that cos 9 — cos 9, which are the roots of the equation 

g{E,k"')^g{E^Uo,k"'), (28) 

(see the definitions ([H])). Particularly, Eg. ipS)) has no solution for a function g{E, fc™) that is monotonic with respect 
to E, except when Uq — 0, the usual condition of the unit transmission. 

For the case of one- mode propagation in semimetal GRs (fc™ — 2tt/2>) the relation (|23p reduces to the simple linear 
dependence 

5(i?,fc") = -i?/2, (29) 

which corresponds to the absence of perfect transmission. In cases of one-mode propagation in semiconducting GRs 
(/c™ ^ 27r/3), due to the oddness oi g with respect to E, we confine our attention to the properties of g{E, fc™) inside 
the positive-energy region E > Q. Differentiating g{E^ k™) with respect to E, 

dgiE.k"') 1 1 - 4cos2 !^ , , , 

= — + — ^ = 0, (fq dE^ > 0, 

dE 2 2E^ ' ^' 

one can compute the point where the function attains its minimum value 



E^E^,nik"') = ^l-icos^—, (30) 

only for transverse wave numbers in the region 27r/3 < A:™ < tt. For the case when < k™ < 2tt/3, we obtain 
dg/dE ^ and d^g/dE'^ < 0. Hence g{E,k"^) is convex and monotone decreasing. Figure [7] presents the in-band 
behaviour of g{E, fc™) with respect to E for the above mentioned fc™-domains: fc™ = 27r/3 — 0.1(< 27r/3), fc™ ~ 2tt/3, 
and A;™ = 27r/3 -I- 0.1(> 27r/3). 

As follows from the relation ([28|) . the transmission probabilities (|2T|) and (|27p reach unity at the points where the 
curve g{E, A:'") crosses g{E—UQ, A;™). Such intersections are possible even for [/q <C 1 in the region 27r/3 < A;™ < tt (see 
Figure [7]), where the resonance occurs near the energy level (|30p . In the case when < k™ < 2tt/3, unit propagation 
takes place if the value of Uq exceeds the energy gap A. To show this, we plug in the relation into Eq. and 
get 

A;™ 

E{E-Uo)^x, x = l-4cos2— . (31) 
The unit transmission points are real- valued roots of the Eq. pip . 

E'^-' = Uo/2 ± V((7o/2)2+x. (32) 
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FIG. 8: Resonance energies iS"^"^^ (solid lines) versus Uq for transverse wave numbers fc™ = 2tt/3 ± 0.1. Dashed horizontal lines 
correspond to the boundaries of the positive- and negative energy subbands for the left region, where A is the energy gap for 
given fc". 



In the interval < /c™ < 2tt/3, x is negative, and the transmission resonances exist only if \Uq\ exceeds 2|x|^/^. 
Expressing x in terms of the energy gap A (see Eg. psp ). we obtain the following inequality for transmission resonances 
in this region 



|f/o| > 4v/A(A + l). 



(33) 



The energies of the transmission resonances versus J7o are plotted in Figure [51 If 27r/3 < fc™ < tt, the indicated 
region |J7o| ^ 1 contains two resonances, one of which belongs to the valence band, and another to the conduction 
band. In the case when < /c™ < 27r/3, two unit transmissions are observed if the condition ([33]) holds. The resonance 
positions depend on the sign of C/o- Both of them belong to the valence band if J/q > and to the conduction band 
if Uo < 0. 

The case of one- mode propagation in semimetal ribbons (fc™ = 27r/3) deserves additional attention. The linear 
dependence of g{E,k™') with respect to E (given by Eq. (|29p ). makes unit propagation impossible. This shows that 
the Klein parado3*i^ does not hold in armchair GRs. Specifically, plugging in ([^ and ([^ into (PT|) produces 



T = 



which can be reduced to 



(^V^^~W+^i-{E-Uory + 



T « 1 - 



(34) 



(35) 



for£; = Oand |;7o| < 1. 

It is worth noting that the through-step expression p4p coincides with its analog for electron transmission through 
a potential step in an undimerized linear chain of identical atoms (see Appendix B). Specifically, the through-step 
coefficient (1341) in the vicinity of the neutrality point i? = is exactly the same as the transmission coefficient through 
a step potential formed in the linear chain when the electron energy is located near the corresponding band center. 



X. LOW-ENERGY LIMIT 



Most theoretical studies on the electronic properties of all-carbon honeycomb lattices are based on the k • p ap- 
proximation in the vicinity of the zero-energy point, leading to the two-dimensional analog of the relativistic Dirac 
equationii^ii^iii^. Rewriting the energy dispersion ^ 



E' 



1 ± 2 cos ■ 



=F 4 cos — 1 — cos — 
2 V 2 



we see that in the low-energy limit only the branch 



= ( 1 - 2 cos — j + 4 cos — f 1 - cos — 



(36) 
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of the dispersion relation takes part in electron transport. In terms of the deviation of k = {k",k™) from the 
zero-energy point (0, 27r/3), 

= 3A:,ac-C, A:" = 2^/3 + VSkyac^c, 
we can expand (|36p into series \kx{y)\ac-c ^ 1 up to and including the quadratic terms. Then 

E{k,,ky)^±^^^kl + kl (37) 



which is similar to the cone-like form of electron dispersioni^^ii^. In the same way, the formula ^ for inter-cell 
electron phase shift 6 reduces to 

,« l-2cosVe^'= zsgn(i?)A±^ + 0(fc,,fc^) + ..., (38) 



where 



E /i.2 4_ 1,2 



0{kx,ky) . sgn(i?) + '^'-'^ + (39) 



Retaining only the first term in Eq. ([38)) . we come to 



e'" = i sgn{E) e"^, = tan"! ^, (40) 

where is the known phase factor from low-energy theoryi^ii^. Expressing the transmission coefhcients ([2T|) and ([27|) 
in terms of 0, (/) = — > fc^;), we obtain the following approximate expressions 

j^stcp ^ 4|COS0 COS 01 ^^^^ 

(I COS 01 -I- I COS 01) -I- [sgn(i? — Uq) sin0 — sgn(i?) sin0] 



)COS 



2 



Tbar ^ ^ ^ (42) 

cos^ 0cos^ -I- [sgn(i? — Uo) sin0 — sgn(i?) sin0] sin^(fcj:£>) 

for step-like and barrier-like potentials in the low-energy limit. The length D ~ 3ac-cN in (|42ll is the width of the 
barrier, as depicted in Figure [51 

The through-barrier coefficient coincidents with the expression^- obtained in the framework of the Dirac 
approach (see details in Appendix [C)) . This coincidence is not surprising. However, no "Dirac" analog of the through- 
step coefBcient (PT|) has been obtained so far. 

Eqs. (|¥T|) and (|^^ show that step-like and barrier-like potentials become transparent in the case of the so-called 
"normal incidence", ky = (or = = 0), which produces the Klein paradox in graphenei^. However, such 
transparency of potentials does not occur if exact tight-binding solutions are written in the low-energy limit. Indeed, 
the expression for cos 9 is determined by the general expansion ((38|l , 



cos 9 = - - 



U2 I jU2 



ky + ^{3kl + kl)+...]. (43) 



At fcj^ = 0, the first term in (|43p vanishes and the expression for cos 9 is described by the next non- vanishing term 
0(kx, ky), 

cos0 = -^sgn(i;)fc,ac-c--|, |S| < 1. (44) 

Inclusion of the next expansion terms in the expressions for cos 9 and cos 9 resembles extending the standard k • p 
method by including higher order k • p terms (see Reffl^ and Eqs. (2.2) and (4.1) therein). For instance, the effect 
of small backscattering in graphene structures was predicted for the first time within the framework of the extended 
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FIG. 9: Through-step transmission coefficients versus the applied potential energy Uo for three different transverse wave 
numbers k"": a) fc" = 27r/3 - 0.018 and Ef = 0.02; b) k"" = 2n/3, Ep = 0; c) = 27r/3 + 0.018, Ef = 0.02. Solid lines 
correspond to the exact solution (|21|) . dashed lines to the approximate dependence (|41|) . The inserts show the same curves in 
the low- voltage region |i7o| < 0.1. Two arrows point to the peak and dip of the exact transmission coefficient (|21|) . 

model^^. Thus, we conclude that the relativistic approach based on the standard k • p scheme should not be applied 
without amendments for describing electron transport through the gapless modes. The next non-vanishing term of 
the expansion ([^5]) "disables" the Klein paradox: the expression {cos 9 — cos 9)^ entering Eqs. (PT|) and ([TT]) becomes 
quadratic on Uq, which results in the same order deviation of the exact transmission coefficients (PT|) and (E7)) from 
unity (see Eq. ((55)) ). 

Evidently, since \ky\ increases proportionally to \k™ — 27r/3|, we can restrict attention only to the first term in (|43|) 
if \ky\ is large enough, which makes the analytical solutions (|2ip and (|27p similar to the approximate dependencies 
(HTjl and for the Dirac formalism. 

XI. MODELING AND DISCUSSION 

Now we compare the exact tight-binding solutions ([2T|) and ([27| with the approximate expressions ([4T|) and (|42p . 
Figure [S] represents ?7o-dependence of the through-step transmission coefficients (l?T|) and (|¥T|l for three different values 
of fc™: (a) fc™ = 27r/3 - 0.018, (b) fc™ = 27r/3, and (c) fc"* = 27r/3 + 0.018. We assume the electron energy E is fixed: 
(a,c) E — Ep = 0.02 and (b) E = Ep = 0. Simple estimations show that this choice leads to single-mode propagation 
in graphene ribbons having (a) M = 58, (b) M = 59, and (c) A/" = 60 elementary cells in the transverse direction. 
These curves are very closed to each other when |[/o| < 0.1, as depicted in Fig. O Thus, the approximation (|¥T|) for 
the low- voltage region is numerically indistinguishable from the exact solution (pij) . 

When \Uq\ < 1, the exact tight-binding solution (PT|) for fc™ — 2tt/3 + 0.018, unlike the approximate solution (|4ip . 
reveals two extreme points, as indicated in the Figure [9j;. The existence of such transmission peculiarities follows 
from the definition an the following chain of identities 

dT _dT d9 _ dT dcos9 _ dT dg{E ~ Uq) _ dT dg{E) 

dUo d9 dUo dcos9 dUo dcos9 dUo dcos9 dE e^e~Uo 

Thus, the additional peak and dip of the transmission coefficient are the extreme points of g{E,k"^ > 2tt/3), 
±y/l - 4cos2(fc™/2), as discussed in Sec. HXl 

It is also important that the exact dependence ((2T|) for k™ ^ 27r/3 predicts unit transmission at applied electrostatic 
energy J7q^* = Ep — x/Ep, which is the solution to pTjl with respect to Uq. However, in the case under consideration 
these resonances are out of the indicated region: Uq'^'^ > 1 for k™ < 27r/3, and Uq'^'^ < —1 if A:'" > 27r/3. 

Figure [TOl depicts the through-barrier transmission coefficients (|27)) and (|42l) for the same values of A;™ and Ep, as 
Fig. m To make the curves more distinguishable, we consider a narrow barrier, N — 10. One can see the series of 
pronounced peaks associated with the resonance conditions k"'N = kxD = nl (1 = 0,1,.. .). A visual comparison of 
the approximation (021) with the exact solution ([77]) shows that they are closed to each other in the low-voltage region. 
However, there are some discrepancies in the resonance positions, especially in the depth of dips outside the low- voltage 
region. Discrepancies of this kind are also observed between the exact and approximate through-step solutions (j2ip 
and (|4ip over the entire C/g-region. Formally, if the barrier-induced system is considered as a combination of two 
scatterers - interfaces between the external and disturbed leads, the through-barrier probability T^^'^ can be written 
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FIG. 10: Through-barrier transmission coefficients (|27p (solid line) and (|42p (dashed line) versus Uq calculated for a) fc™ = 
27r/3 - 0.018 and Ep = 0.02; b) fc"* = 2tt/3, Ef = 0; c) fc™ = 2n/3 + 0.018, Ef = 0.02. The grey rectangle is the region of 
tunneling through a gap which is not studied here. 



down in terms of the through-step transmission T^*-°p, The probability T^^'^ varies within the interval defined by (32) 
from Refi^ 

< T^''' < 1, 

which can also be observed from the behaviour of the corresponding through-step coefficient depicted in Fig. [S] 

Obviously, as |fc™ — 27r/3| increases the exact (tight-binding) and approximate solutions for through-step and 
through-barrier transmissions in the low-voltage region (|C/o| <C 1) are closer to each other. This also follows from 
Sec. IXl specifically from the expansion (|43)) . 

Now we consider propagation through a gapless band (/c™ ~ 27r/3). Despite visual coincidence of the exact tight- 
binding solutions with unity in the low- voltage region (Figs. andflOb) the dependences pT|) and ([27|) do not reveal 
unit propagation for the whole low-voltage region. Thus, the Klein paradox, which strongly predicts that T = 1 for 
one-mode propagation through any high and wide potential barriers in semimetal ribbons, fails since the through-step 
and through-barrier coefficients 2^stcp(bar) (jjffgp fj-om unity though the corresponding deviations are indistinguishable 
at experimentally reliable values \Uo\ ^ 1. 



XII. CONCLUSIONS 



In the framework of the tight-binding model, we present a full and closed description of electron transport in 
armchair GRs. Starting with the tight-binding Schrodinger equations for a 2D honeycomb lattice, we identify the 
boundary conditions, wave function and propagation velocity of an electron in armchair GRs. We obtain analytical 
expressions for the through-step and through-barrier transmission coefficients and demonstrate that new type trans- 
mission resonances exist. Such resonances occur when cos 9 = cos 9, where 9 and 9 are the inter-cell electron phase 
shifts in the corresponding regions. We show that the number and positions of the resonances are strongly dependent 
on the electron transverse wave number fc™. In particular, at fc™ = 27r/3 (propagation through a gapless band), the 
resonances are absent, and the through-step transmission coefficient reduces to the relation known for through-step 
propagation in a linear chain of identical atoms. For the low-energy limit, the deviation of the through-step trans- 
mission coefficient from unity is proportional to the square of applied potential energy. A similar deviation is also 
observed for the through-barrier transmission coefficient. 

The discrepancy between our tight-binding result and the Klein paradox in graphene becomes evident as a result of 
expanding the expressions for cos^? and cos^ in the vicinity of the zero-energy point. We show that in the low-energy 
limit these expressions come to negligibly small terms, which are ignored in the well-known k • p method leading 
to the relativistic Dirac equation. These negligibly small terms "destroy" unit propagation and produce a small 
backscattering proportional to Uq. 

The presented analytical results are in complete agreement with the results of the numerical computations made 
in the paper. 
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APPENDIX A: COMPUTATION OF ELECTRON FLUX 



In this appendix, we compute the flux related to electron wave propagation in armchair graphene ribbons. We start 
with the set of non-stationary Schrodinger equations 



'4^^n,m,p — ['0n,m,A ~t~ ^n,m.r ^ (f ^m,l)'0n,m — l,r] ; 



l<m<N, 



Multiplying the equations for ipn,m,a by V'n m a {m = 1, . . . ,Af + 1) and adding the obtained relations to the complex 
conjugated counterparts, and then summing them up over all carbon atoms in unit cell n, we get the following 



where 



dt 



— 1 jright 



(Al) 



jleft ^_ [V'„-l,„,rCm,i - C-l,m,rV'n,m,;] 



^^ 



and 



jright 



n+l,m 



m—1 



The left side of Eq.(Al) is the rate of change of the electron probability in unit cell n. The right side of the equation 
gives the flux entering to the unit cell. Thus, the flux J outgoing from unit cell n to the right direction equals 



J = _ J'''9''* 

This formula is used in the main text. 



i(3 



n+l,m,^ V'n 



{A2) 



APPENDIX B: THROUGH-STEP TRANSMISSION IN A LINEAR CHAIN OF ATOMS 



In this section, we calculate the through-step transmission coefficient for electron moving in a linear chain of identical 
atoms (see Fig. [TT|) . For the sake of convenience, we keep the same notation for E = E/\(3\ and Uq = eVb/|/3|. 
Following Fig. Ill[ we need to determine the solution to the following system of linear Schrodinger equations 



left left left 

^u,-E)i.r=^T+^r, 

;ight right right 



(Bl) 



V'n-l+V'n+l, n>l, 



where 



re 



-ikn 



te 



ikn 



(B2) 
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FIG. 11: Model of the linear chain of atoms coupled by the electron-transfer integral /3. As in the main text, atoms with n < 
correspond to the left lead, and n > 1 - to the right lead. 



r and t are the amplitudes of the reflected and transmitted waves, respectively, k and k are the non-dimensional wave 
numbers in the left and right regions. One can see that the first and forth equations in (Bl) are satisfied identically if 

cosfc = — — , cosfc = 2~^' i^'^)- 

Plugging in (B2) into the second and third equations from (Bl) gives the system of two linear equations 

1 + r = t, e + e r — e t. 

The solution to the system is 

gifc — e''^ 2i sin k 

t 



Thus, the transmission coefficient determined as T = |i| sinfc/sinfc can be represented as 

^ 4 sin k sin k , , 

T= = = 2 2- (B4) 

4 sin k sin fc -h (sin k — sin k) + (cos k — cos k) 

Since the values of sin A: and sin A: are positive in the corresponding bands (0 < (fc,fc) < tt) one can come to an 
expression similar to (|34p for the transmission coefhcient. 

Therefore, independently on the Fermi position inside the band, electron propagation through a potential step in 
a graphene ribbon is described by the same through-step relation as in an undimerized ID chain with the same ratio 
Ua = eVom. 



APPENDIX C: RELATIVISTIC DIRAC SOLUTION FOR THE ELECTRON REFLECTION 
PROBABILITY THROUGH A POTENTIAL BARRIER. 



The reflection amplitude from a potential barrier obtained through the Dirac formalism is described by Eq. (3) in 
Ref-i^ 

2ie'^'^ sm(kxD) [sin (j) ss' sin (f)] (CI) 
ss' [e"*'^^-^ cos{(f> + (j>) + e^^'^-^ cos(0 — i^)] — 2isin{kxD) ' 

where s = sgn(E) and s' — sgn(£' — Uq). Equivalently, 

ie'*^ sin(fc2;_D) [sin </> — ss' sin (j)] 
ss' cos{kxD) cos cos <j) + i sin{kxD)[ss' sin(/)sin (j> — 1] 

Since 

(1 — ss' sin sin (j))'^ = (sin (j) — ss' sin (f>)'^ + cos^ (p cos^ (j), 

the reflection probability R = |rp is given by 

^ ain^ (kxD) [sin <j) ~ ss' sin (j)]"^ 

cos^ 0cos^ (f> + sin^(fcj:-D)[sin (j) — ss' sin0]^ 
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Thus, the transmission probabiHty T = 1 — R through the potential barrier is exactly the same as our approximate 
relation (li^ . 
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